Preprint typeset using I4TgX style emulateapj v. 6/22/04 
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ABSTRACT 

There is growing evidence that two classes of high-energy sources, the Soft Gamma Repeaters and the 
Anomalous X-ray Pulsars contain slowly spinning "magnetars", i.e. neutron stars whose emission is powered 
by the release of energy from their extremely strong magnetic fields (> 10'^ G). We show here that the enor- 
mous energy liberated in the 2004 December 27 giant flare from SGR 1806-20 5 x lO**^ erg), together with 
the likely recurrence time of such events, requires an internal field strength of ^ 10"' G. Toroidal magnetic 
fields of this strength are within an order of magnitude of the maximum fields that can be generated in the core 
of differentially-rotating neutron stars immediately after their formation, if their initial spin period is of a few 
milliseconds. A substantial deformation of the neutron star is induced by these magnetic fields and, provided 
the deformation axis is offset from the spin axis, a newborn fast-spinning magnetar would radiate for a few 
weeks a strong gravitational wave signal the frequency of which 0.5-2 kHz range) decreases in time. The 
signal from a newborn magnetar with internal field > 10'^^ G could be detected with Advanced LIGO-class 
detectors up to the distance of the Virgo cluster (characteristic amplitude h^- ^ 10"^'). Magnetars are expected 
to form in Virgo at a rate ^ 1 yr"' . If a fraction of these have sufficiently high internal magnetic field, then 
newborn magnetars constitute a promising new class of gravitational wave emitters. 

Subject headings: gravitational waves — stars: magnetic fields — stars: neutron — stars: individual(SGR 
1806-20) 



1. INTRODUCTION 

The two classes of pulsating high energy sources that 
are believed to host magnetars, the Soft Gamma Repeaters, 
SGRs, and the Anomalous X- ray Pulsars, AXPs, share a 
number of propertie s (see e.g. Mereehetti & Stella ( 1995|; 
IKouveliotou et all il99M; Wood s & Thomp son (2004)). 
They have fairly long spin periods in the ~ 5 10 s range, 
spin-down secularly with timescales ~ 10"* 10^ yr, do not 
have a companion star and are in some cases associated to 
supernova remnants with ages ~ 10^ 10^ yr Rotational 
energy losses are 1^2 orders of magnitude too low to explain 
the persistent emission of these sources, which is typically 
~ 10^"*^ 10^^ erg/s. Both AXPs and SGRs show periods of 
activity in which recurrent short bursts (<?; 1 s) are emitted, 
with peak luminosities in the ^ 10^^ - 10'*' erg/s range. 
The initial spikes of giant flares are second-long events of 
exceptionally high luminosity, 3 to 6 orders of magnitude 
larger than that of the brighest recurrent bursts. Giant flares 
are rare, only three have been observed so far, and might 
well be specific to SGRs. Given the highly super-Eddington 
luminosities of recurrent bursts and, especially, giant flares, 
accretion power is not an option. 

The magnetar model envisages that SGRs and AXPs 
are powered by the release of energy from their ex- 
tremely high magnetic fields (Duncan & Thompson' Il992t 
[Thompson & Duncan 1993, 1995, 1996, 2001). The model 
has been largely successful in interpreting the unique fea- 
tures of these sources. According to it, a twisted, mainly 
toroidal magnetic field characterizes the neutron star inte- 
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rior (B > 10'^ G), only partially threading the crust. The 
emerged (mainly poloidal) field makes up the neutron star 
magnetosphere, with dipole strenghts of ~ few x 10'^ G, 
as required to generate spi n -down at the observed rat e 
dThompson & DuncanI Il993t iThompson & Murravl l2001h . 
Energy is fed to the neutron star magnetosphere through 
Alfven waves driven by local "crustquakes" of various am- 
plitude and producing recurrent bursts. Giant flares may 
result from large-scale rearrangements of the core mag- 
netic field or catastrophic instabilities in the magnetosphere 
(Thompson & Duncan 2001; Lvutikov 2003). These events 
occur more rarely and lead to the sudden release of very large 
amounts of energy. Most of this energy breaks out of the mag- 
netosphere in the form of a fireball of plasma expanding at 
relativistic speeds and giving rise to the initial spike of gi- 
ant flares. The oscillating tail that follows this spike, display- 
ing up to ~ 50 cycles of the neutron star spin, is interpreted 
as due to that part of the flare energy that remains trapped 
in the magnetosphere. The total energy released in this tail 
(~ 10'*"' erg in all three events detected so far) provides an 
independe nt estimate of ^ lO'"' G for the e xternal field of 
magnetars ( IThompson & DuncanI 1 9951 l200l . 

2. THE INTERNAL B-HELD OF MAGNETARS 

The 2004 December 27 giant flare from SGRl 806-20 pro- 
vides a new estimate of the internal field of magnetars. About 
5 X 10"*^ e rg were released in the ~ 0.6 s long mai n spike of 
this event dTerasawa et al.l20()^lHurlev et alJ20()5l) i.e. a two 
decade higher energy than that of the other giant flares ob- 
served so far, the 1979 March 5 event from SGR 0526-66 
(|Mazets et al. 1979) and the 1998 August 2 7 even t from SGR 
1900H-14 (Hurlev et al. 1999: Feroci et aniT999l) . Only one 
such powerful flare has been recorded in about 30 yr of mon- 
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itoring from the ~ 5 known magnetars in SGRs. The recur- 
rence time in a single magnetar implied by this event is thus 
about 150 yr (if giant flares were emitted also by AXPs, 
this recurrence time would increase by a factor of ^ 2, as 
the number of AXPs and SGRs in our galaxy is compara- 
ble). The realisation that powerful giant flares could be ob- 
served from distances of tens of Mpc (and thus might repre- 
sent a sizeable fraction of the short Gamma Ray Burst popu- 
lation) motivated searche s for 2004 Dec 27-Uke events in the 
BATSE GRB database ( Lazzati et al. 2005; Po pov & Steriil 
|2005). The upper limits on the recurrence time of powerful 
giant flares obtained in these studies range from r ^ 130 to 
600 yr per galaxy, i.e. ^ 4 to 20 times longer than the time 
inferred from the 2004 Dec 27 event. A detailed discussion 
of the merits of these estimates is beyond the scope of this 
Letter However, we note that a 2004 Dec 27 -like event in the 
Galaxy could not be missed, whereas several systematic ef- 
fects can reduce the chances of detection from large distances 
(see e .g. the discussion in Lazzati et al. (2005); Nakar et al. 
(|2005)). Rather than regarding the 2004 Dec 27 event as sta- 
tistically unlikely, one can evaluate the chances of a recur- 
rence time of hundreds of years, given the occurrence of the 
2004 Dec 27 hyperflare. We estimate that, having observed a 
powerful giant flare in our galaxy in ^ 30 yr of observations, 
the Bayesian probability that the galactic recurrence time is 
T > 600 yr is ^ 10"^, whereas the 90% confidence upper limit 
is r ^ 60 yr. These valeus are derived for a uniform prior' on 
the rate 1 /t. We thus favor smaller values and assume in the 
following T ^ 30 yr. 

In a timescale of '--^ 10'* yr (that we adopt for the SGR life- 
time) about 70 very powerful giant flares should be emitted 
which release a total energy of 4 x 10'*^ erg. We note that if 
the powerful giant flares' emission were beamed in a fraction 
b of the sky (and thus the energy released in individual flares 
a factor of b lower), the recurrence time would be a factor of 
b shorter. Therefore the total energy released by giant flares 
would remain the same. For this energy to originate in the 
magnetar's B-field, this must be > 10'^^ G. This value should 
be regarded as a lower limit on the initial internal field of a 
magnetar. Firstly, not all the the magnetic energy is released 
through powerful giant flares. The magnetar model predicts 
a conspicuous neutrino luminosity, resulting from ambipolar 

diffusion, of L^, ~ 4 x 10^^ t^'^^'' erg s"', where t4 is the mag- 
netar age in units of lO'* yr (Thompson & Duncan 1996). Inte- 
grating over f4 = (0.01 ^ 1) yields a neutrino energy output of 

4 X 10'**' erg. Including this energy gives an internal mag- 
netic field of > 10'^^ G. Secondly, ambipolar diffusion and 
magnetic dissipation are expected to take place at a faster rate 
for higher values of the field (Thompson & Duncan 1996). 
Therefore estimates of the B-field based on present day prop- 
erties of SGRs (age 10"* yr) likely underestimate the value 
of their initial magnetic field. 

Very strong internal B -fields are expected to be generated 
in a differentially rotating fast spinning neutron star, subject 
to vigorous neutrino -driven convection instants after its for- 
mation cDuncan & Thompson 1992,) . If the initial spin period 
Pj is in ~ 1 2 ms range (comparable to the overturn time 
of convective cells), an efficient a - dynamo is powered by 
fluid helycal motions, which amplifies the field in a predom- 
inantly toroidal configuration with a coherence length com- 

' For a uniform prior on log( 1 /r) and a (conservative) lower cutoff of 
1/(10V). 

we obtain t > 600 yr with probability 0.04 and a 90% confidence 
upper limit of r ~ 270 yr. 



parable to the star size. Differential rotation provides a very 
large free energy reservoir, ^ 10^^ (P,/ 1 ms)"^ erg, that can be 
converted into a magnetic field of up to ^ 3 x 10''(P, / 1 ms)"' 
G ( Duncan & Thompson 1992). However the efficiency of 
the amplification is likely limited by the dynamical influence 
(backreaction) of the field which opposes convective motions. 
This limitation could be circumvented, however, if smooth- 
ing of angular velocity gradients by magnetic stresses is very 
efficient on a timescale shorter than the initial neutrino cool- 
ing time (which sets the timescale for convective instability 
to occur); in this case, a field of several x 10'^ G can be gen- 
erated in ma gnetars that are born with spin periods of a few 
milliseconds fThompson & Duncan ll993HDunca nll99a) . 

In summary, unless powerful giant flares such as the 2004 
December 27 event are much rarer than the rate implied by 
having already detected one of them, magnetars must possess 
internal fields of lO'^ G or higher (values up to ^ 10'^ G 
cannot be ruled out). Based on current models, fields of this 
order can be generated only if magnetars are born with spin 
periods of a few ms. In the following we explore the conse- 
quences of this for the generation of gravitational waves from 
newborn magnetars. We parametrize their (internal) toroidal 
field with B,,i6.3 = B,/2 x 10'^ G, (external) dipole field with 
Bd.i4 = Bd/ 10^* G and initial spin period with P, 2 = Pi/ (2 ms). 

3. MAGNETICALLY-INDUCED DISTORTION AND GRAVITATIONAL 
WAVE EMISSION 

The possibility that fast-rotating, magnetically-distorted 
neutron stars are conspicuous sources of gravitational radi- 
ation has been discussed by several authors (see e.g. the re- 
view by Bonazzola & Marck ( 1994)). The path towards a bet- 
ter understanding of the early evolution and gravitational ra- 
diation emission_groperties of such neutron stars was layed 
out by ICutled (|2003), who discussed also the possibility of 
detecting the gravitational radiation signal from a newborn, 
^ 10'^ G neutron star in our Galaxy. If this value for the in- 
ternal field pertains to ordinary neutron stars, then such events 
might be expected at the galactic rate of neutron star forma- 
tion, i.e. once every several tens to hundreds of years (see e.g. 
'Arzoumani an et akl J2002h ). The formation rate of magnetars 
is substantially lower than that (see Section 4), making the 
chances of witnessing the birth of a galactic magnetar in the 
near future very slim. In the following we show that for the 
range of magnetic fields discussed in Section 2, newly born, 
millisecond spinning magnetars are conspicuous sources of 
gravitational radiation that might be detectable up to the dis- 
tances of the Virgo cluster 

The anisotropic pressure from the toroidal B-field deforms 
a magnetar into a prolate shape, with ellipticity ~ -6.4 x 
10~^(< B^i6 3 >), where the brackets indicate a volume- 
average over the entire core (Cutler 2002). In general, the 
symmetry axis of the magnetically-induced deformation will 
not be coaligned with the rotation axis. Hence, in the neu- 
tron star frame the angular velocity vector undergoes free 
precession with period Pprec — P/^b around the (magnetic) 
symmetry axis, where P is the spin period (Jones 2OO4 
iJones & And ersson 2002 ) (the rotationally-induced distortion 
is not relevant here, because it is always aligned with the in- 
stantaneous spin axis). As long as the axis of the magnetic 
distortion is not aligned with the spin axis, the star's rotation 
will cause a periodic variation of the mass quadrupole mo- 
ment, in turn resulting in the emission of gravitational waves, 
GWs, at twice the spin frequency of the star. 

A crucial issue is whether the misalignement of the mag- 
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netic axis can be maintained in spite of the damping of free 
precession caused by the star's internal viscosity. In a prolate 
ellipsoid the rotational energy is minimized when the moment 
of inertia along the axis of t he "frozen-in" distortion, 7^, is 
orthogonal to the spin axis (J ones & And ersson 2002; Cutler 
|20O2). Hence, viscous damping of precession drives a pro- 
late neutron star towards a geometry that maximizes the time- 
varying mass quadrupole moment and, thus, GW emission. 

As the power emitted in GWs scales as oc P~^, the GW sig- 
nal, for a given internal B-field, depends critically on the ini- 
tial spin period and its variation thereafter. The spin evolution 
of a newborn magnetar is determined by angular momentum 
losses from GWs, electromagnetic dipole radiation and rela- 
tivistic winds. According to Thompson et al. (2004), the lat- 
ter mechanism depends critically on the value of the external 
dipole and is negligible except for fields > (6 7) x 10'^ G, 
a regime that we do not discuss here. We consider in the fol- 
lowing the spin evolution of a newborn magnetar under the 
combined effects of GW and electromagnetic dipole radia- 
tion. When both mechanisms are taken into account, the mag- 
netar angular velocity lo = 2-k/P evolves according to 



LO = - 



GW 5 



(1) 



where = {B]R^) / {3Ic^) and Kcw = (128 /5)iG/c^)Iel, with 
R the neutron star radius, G the gravitational constant and c 
the speed of light. From this equation, we derive a spin-down 
timescale of 



Tsd = L0/(2lj) ■ 



10 P^(Bi,4+l-15Bf,i6.3^^')" d 



(2) 



According to mode ls of ne utron star internal viscosity 
(|Alpar & Pines 1985: lAlpar & Sauls 1988), orthogonaHzation 
of the magnetic axis is expected to take place in n ~ 10^ pre- 
cession cycles (Cutler 2002). This translates into a timescale 
of Tort = nPi/es ^ 0.4Bj"[g 3P, 2 d. The condition for the newly 
formed magnetar to become an orthogonal rotator before loos- 
ing most of its spin energy is 



Tsd 
Tort 
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(3) 



If condition (|3} is met, the magnetar quickly becomes a max- 
imally efficient GW emitter, while its spin period is still close 
to the initial one. The strain is 



(4) 



where the distance d2Q = d/(2Q Mpc) is in units of the Virgo 
Cl uster distance and the an gle-averaged strain is that given 
bv llJshomirskv et alJ (l20()0ll . We estimate the characteristic 
amplitude, he = /lA^'/^, where ~ Tsd /Pi is the number of cy- 
cles over which the signal is observed. Using equation (|4} we 
obtain 
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Under the conditions discussed above, strong GW losses are 
not quenched immediately after the magnetar birth but rather 
extend in time, typically from days to a few weeks, before 
fading away as a result of the star spin-down. 

In order to assess the detectability of this GW signal we 
compute the optimal (matched-filter) signal-to-noise ratio 



= 2 



\h(f)? 
If, S,(f) 



df 
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by extending the formulae bv lOwen &Lindbloml ( l200l and 
Cutler (2OO3I1 to include both the electromagnetic and GW 
torques. In the previous equation fi = 2/f, w 2 kHz and 
ff = 2/Pf ^ 500 Hz are the initial and final values of the fre- 
quency of the signal. h(f) is the Fourier transform of the GW 
strain at the detector (Eq. (0 divided by 5'/^ in order to ac- 
count for averaging over source sky position) and Shif) the 
one-sided noise spectral density of the instrument. In the rel- 
evant frequency band 500 Hz - 2 kHz we approximate Shif) 
according to the current baseline performance of Advanced 
L1GO:5a(/) = 2.1 x 10-'*^(//lkHz)2Hz"' (Thorn?2000). We 
com pute hi f) using the station ary phase approximation (see 
e.g. lOwen & LindblomI (|2003)) with angle-averaged ampli- 
tude and frequency derivative / = w/tt given by Eqs. (0} and 
([!), respectively. Eq.(|6|l yields the following estimate of the 
optimal signal-to-noise ratio: 



3*^20^?, 16. 3^rf'l4 



In 



+ 21n 



(7) 

where = 2Kd / iii^Kcw)- The limit in which the electromag- 
netic torque dominates corresponds to ^ (see Eq. ([0). 
In this case, the first logarithmic term in Eq. (Q is negligible 
and we recover the result by Cutler ( 2002). Fig.^shows lines 
of constant S/N for a source at c/20 = 1 and selected values of 
the initial rotation period (f, = 1.2,2 and 2.5 ms) in the iB,,Bd) 
plane: GWs from newborn magnetars can produce S/N > 8 
for B, > lO'*"^ G and Bd < lO''*^ G. This is the region of pa- 
rameter space that offers the best prospects for detection as 
we now discuss. 

The optimal (in the maximum likelihood sense) detec- 
tion statistic for long- lived signals is the so-called J^-statistic 
dJaranowski et alJl99 8): its computation involves the correla- 
tion of the data stream with a discrete set of filters that probe 
the relevant space of unknown parameters: the sky position 
(2 parameters) and 2 parameters that control the evolution 
of the GW phase (see Eq.(n]i). A very conservative upper- 
limit on the required number of filters involved in a search 
can be evaluated ig noring correlatio ns among the parame- 
ters: 100 filters dBradv et alJ 119981) are needed to cover 
the ~ 40deg^ area of the Virgo cluster for a coherent inte- 
gration time T - 10*' s, and - iT ~ 10'^ filters for the 
phase parameters, giving a total of ^ 10^^ templates. For this 
number of trials, we estimate that a signal with S/N ^ 12 can 
be detected with 1% false alarm and 10% false dismissal (if 
all the parameters of the signal were known the correspond- 
ing S/N would be « 4.6). Of course the actual sensitivity 
limit of such a search is likely set by the available computa- 
tional resources. In practice, the search for GWs from new- 
born magnetars is probably best carried out using a "hierar- 
chical search strategy" where coherent and incoherent stages 
are alternated in order to achieve (quasi-)optimal sens itivity 
at affordable computational costs (iBradv & Creighton 2000; 
Krishnan et al. 2004; Cutler etal. 2005;' Abbott et al. 2005). 
Furthermore, the presence of a trigger based on the detection 
of the newborn magnetar by other means (e.g. the correspond- 
ing supernova) could further reduce the region of parameter 
space that needs to be searched. In addition the sensitivity 
will increase with the operation of the network of GW detec- 
tors that is coming on line. In any case the search for this new 
class of signals represents a challenge that we are investigat- 
ing in greater depth. 
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4. DISCUSSION 

We have shown that the energy liberated in the 2004 De- 
cember 27 flare from SGR 1806-20, together with the likely 
recurrence rate of these events, points to a magnetar internal 
field strength of ^ 10'^ G or greater Such a field is suffi- 
cient to induce a sizeable deformation of the star and is most 
likely generated in a differentially rotating, millisecond spin- 
ning magnetar instants after gravitational collapse. Magnetars 
with these characteristics are thus expected to be very power- 
ful sources of gravitational radiation during the first weeks of 
their life. Prospects for revealing their GW signal depend crit- 
ically on the birth rate of these objects. The three confirmed 
associations between an AXP and a supernova remnant (ages 
in the 10^ ^ 10"* yr range) impli es a magnetar birth r ate of 
> 0.5 X 10"^ yr"' in the Galaxy (iGaensler et al.llT999l) . The 



chances of revealing the signal from a newborn magnetar in 
our Galaxy are thus very low. A rich cluster like Virgo, con- 
taining 2000 galaxies, is expected to give birth to magne- 
tars at a rate of ^ 1 yr"' . A fraction of these magnetars might 
have sufficiently high toroidal fields that a detectable GW is 
produced. A slowly evolving periodic GW signal at ~ 1 kHz, 
whose frequency halves over weeks, would unambiguously 
reveal the early days of fast spinning magnetars. We conclude 
that newborn, fast spinning magnetars provide a class of GW 
emitters over Virgo scale distances that might well be within 
reach for the forthcoming generation of GW detectors. 
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log B,(G) 

Fig. 1 . — Lines of constant S/N for selected values of the initial spin period in the internal toroidal magnetic field, B,, and external dipole field, Bi, plane for a 
source at the distance of the Virgo Cluster (1^20 = !)■ Solid, dashed and dotted curves correspond to an initial spin period of P, =1.2, 2 and 2.5 ms, respectively. The 
calculations take into account the time required for the toroidal magnetic field axis to become orthogonal to the spin axis. Note that according to Thompson et a^ 
120041) . strong angular momentum losses by relativistic winds set in and dominate the spin down for Bj > 6-^1 X lO'"* G; the curves for such values of should 
thus be treated with caution. 



